Line shape of the fiH(3p — Is) hyperfine transitions 
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The (3p — Is) X-ray transition to the muonic hydrogen ground state was measured with a high- 
resolution crystal spectrometer. A Doppler effect broadening of the X-ray line was established which 
could be attributed to different Coulomb de-excitation steps preceding the measured transition. The 
assumption of a statistical population of the hyperfine levels of the muonic hydrogen ground state 
was directly confirmed by the experiment and measured values for the hyperfine splitting can be 
reported. The results allow a decisive test of advanced cascade model calculations and establish a 
method to extract fundamental strong-interaction parameters from pionic hydrogen experiments. 



PACS numbers: 36.10.-k 

A series of experiments has been conducted at the 
Paul Scherrer Institut (PSI), Switzerland, to extract the 
isospin separated pion-nucleon scattering lengths from 
the observation of X-ray transitions feeding the ground 
state of pionic hydrogen [J 0, Q • With X-ray energies of 
about 2-3 keV and values for the strong- interaction shift 
and broadening being of about 7eV and leV, respec- 
tively, the use of a high resolution Bragg spectrometer 
was mandatory to reach the envisaged precision on the 
per cent level. 

The experimental difficulties are considerable, espe- 
cially for the isovector scattering length, which is equiva- 
lent to the determination of the strong interaction broad- 
ening. It requires to extract from the measured line 
shape a Lorentzian profile representing the natural width, 
which is convoluted with the spectrometer response and 
several contributions owing to the Doppler broadenings 
from different high velocity states of the exotic atom. 

High velocity states develop in exotic hydrogen atoms 
during the atomic de-excitation cascade. As the system 
is electrically neutral, it may dive deeply into the elec- 
tron cloud of a neighbouring hydrogen molecule. Such 
close collisions strongly influence the atomic cascade. 
The most important of these processes are the so-called 
Coulomb transitions , collision-induced radiation- 

less de-excitations, where the released energy is shared 
between the exotic atom and a normal hydrogen atom 
as recoil partner. During the cascade, acceleration due 
to Coulomb transitions and deceleration by elastic and 
inelastic collisions compete, which results in a complex 
and level dependent kinetic energy distribution. 

Historically, the first evidence for high velocity states 



was found in the charge exchange reaction ir~p — > ir°n 
with stopped pions as Doppler broadening of the time- 
of-flight (TOF) of the monoenergetic neutrons Later 
on detailed studies confirmed that Coulomb de-excitation 
substantially affects the kinetic energy of 7rH [8] and y«H 
atoms even at lowest densities. From the TOF spectra, 
several components were identified and attributed to spe- 
cific Coulomb transitions. Hints for the influence of the 
Doppler broadening in X-ray transitions were identified 
in experiments measuring the strong-interaction width of 
the 7rH ground state [2j, |3j . A correction for the cascade- 
induced broadening is therefore indispensible for a proper 
extraction of the hadronic contribution to the X-ray line 
width. However, the high precision information from the 
reaction 7r~p — > Tr°n cannot be directly transfered to 
X-ray studies. Charge exchange occurs from ns states, 
mainly with principal quantum numbers n = 3 — 5 of the 
tt~P system, whereas the initial states for K X-ray emis- 
sion are np levels. Consequently, the preceding cascade 
steps for the two processes are different and the Doppler 
contributions to the X-ray line shape deviate from the 
ones derived from neutron TOF experiments. 

As the intensities of the X-ray transitions strongly 
depend on the hydrogen density, there was a first ap- 
proach to extract information about the cross sections of 
Coulomb de-excitation from intensity studies. Different 
processes like Stark transitions and external Auger ef- 
fect, however, are overwhelming in their importance for 
the X-ray intensity compared to Coulomb de-excitation. 
In consequence, the kinetic energy T^in of the exotic hy- 
drogen atom was used in earlier cascade codes (the so- 
called Standard Cascade Model: SCM) as a fitting pa- 
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rameter with values around T^in = 1 eV, which explained 
the measured intensities with sufficient accuracy [10| . 

Based on the SCM an extended standard cascade 
model (ESCM) was developed. It is a new approach to 
calculate the Doppler contributions to neutron TOF and 
exotic hydrogen X-ray spectra by taking into account the 
competing processes in each de-excitation step and a ki- 
netic energy distribution at the time of X-ray emission is 
provided An example of such a kinetic energy dis- 
tribution for muonic hydrogen in the 3p state is shown 
in Fig.[TJ Monoenergetic lines corresponding to specific 
Coulomb transitions n — > n' are smeared out because of 
the numerous elastic collisions after the Coulomb transi- 
tion and before X-ray emission. The validity of ESCM 
calculations cannot be tested directly in pionic hydrogen 
as the strong-interaction broadening completely masks 
the fine details of the various Doppler contributions. 

Muonic hydrogen as a purely electromagnetic twin sys- 
tem to pionic hydrogen offers itself as an ideal candidate 
for the direct observation of the Coulomb de-excitation. 
Ideally, it could be used as a test of the validity of the 
ESCM by reproducing the line shape of muonic hydro- 
gen X-ray transitions in a fitting routine with the predic- 
tions of the ESCM as input. The theoretical values for 
the ground state hyperfine splitting, and more impor- 
tant, for the relative intensity of the transitions feeding 
the triplet and singlet components should be obtained. 
The ground-state hyperfine splitting is calculated to be 
182.725±0.062meV [12| . A relative statistical population 
of 3:1 is expected for the triplet and singlet components. 

An understanding of the atomic cascade in exotic 
hydrogen is needed in other experimental studies as 
well, e.g., (i) for the precision determination of the 
proton charge radius from the muonic hydrogen 2s-2p 
Lamb shift [3, EH or (ii) in the measurement of the in- 
duced pseudoscalar coupling in muon capture by the pro- 

ton[H[l![3. 

The experiment was performed at the 7rE5 channel 

of the proton accelerator at PSI, which provides a low- 
energy pion beam with intensities of up to a few 10 8 /s. 
Pions of 112MeV/c were injected into the cyclotron trap 
II [l[ and decelerated using a set of degraders optimized 
to the number of muon stops by measuring X-rays from 
muonic helium. Muonic atoms are formed by slow muons 
originating from the decay of almost stopped pions close 
to or in a cylindrical cryogenic target cell of 22 cm length 
and 5 cm in diameter in the center of the trap. The cell 
was filled with hydrogen gas cooled down to 25K at 1 bar 
absolute pressure, which corresponds to a density equiv- 
alent to 12.5 bar at room temperature. X-radiation could 
exit the target cell axially through a 5 /im thick mylar® 
window. 

X-rays from the /iH(3p — Is) transition were measured 
with a Johann-type Bragg spectrometer equipped with a 
spherically bent Si(lll) crystal having a radius of curva- 
ture of 2982. 2±0. 3mm and a free diameter of 90mm [181 . 
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FIG. 1: Kinetic energy distribution of fiR atoms in the 3p 
state for a hydrogen density equivalent to 12.5 bar as pre- 
dicted from the ECSM [ll[ using a) cross sections calculated 
by [llj (dashed curve) and b) recalculated cross sections [l7j 
stimulated by the present experiment (solid curve). Num- 
bers indicate the corresponding Coulomb transitions n — ► n'. 
Note the change of the vertical scale at 10 eV. 



Such a spectrometer is able to measure simultaneously 
an energy interval according to the width of the X-ray 
source when using a correspondingly extended X-ray de- 
tector, in this case a 3 x 2 array of charge- coupled devices 
(CCDs) covering in total 72 mm in height and 48 mm in 
width [19] . A pixel size of 40 /im provides the two dimen- 
sional position resolution necessary in order to measure 
the diffraction image. The background rejection capa- 
bility of CCDs together with a massive concrete shield- 
ing suppresses efficiently events from beam induced reac- 
tions. In total almost 10000 events were collected for the 
/zH(3p-ls) line (Fig.©. 

The spin-averaged fiH('3p — Is) transition energy is cal- 
culated to be 2249.461 ± 0.001 eV with a radiative line 
width of 0.3 ueV [20j. The 3p-level splittings amount 
to a few meV only [211]. Hence, two components with 
identical response functions are sufficient to describe the 
(3p— Is) line. The spectrometer response was determined 
using narrow X-rays from helium-like argon as outlined 
in [22|, [23| . Applying this method to chlorine and sul- 
fur, an extrapolation yields a resolution of 272±3meV 
(FWHM) at the /iH(3p — Is) transition energy, which 
is significantly narrower than the observed line width 
(Fig.©. Details may be found elsewhere 24j . 

As a first trial to include the Doppler broadening 
caused by Coulomb de-excitation the kinetic energy dis- 
tribution given by the ESCM result of [ll[ (Fig.Q]) was 
taken directly as an input for the fitting of the line shape, 
which was done by means of the MINUIT package [25| . A 
comparison to the measured line shape yields a poor re- 
duced x 2 01 Xr = 1-353 only (Fig.O-dashed dotted line). 

It is evident that this ESCM prediction with a weight 
of 36% only for the low-energy component Tkin < 2eV 
underestimates substantially the fraction of X-rays with 
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FIG. 2: Line shape of the ^H(3p — Is) transition as mea- 
sured with a Si (111) crystal in first order. One energy bin 
corresponds to 2 CCD pixels or 37.2 meV. The spectrom- 
eter response (thin solid line) represents the expected line 
shape formed by the two hyperfine components (dashed lines) 
without any Doppler broadening (normalised to the peak 
height after background subtraction). Including Coulomb de- 
excitation by using long standing cross sections [ill ] yields a 
poor description (dashed-dotted) of the line shape. The "best 
fit" in a "model free" approach (see text) in shown by thick 
solid line following the data. A very good fit - being indistin- 
guishable by eye from the "model free" approach - was found 
when using for the ESCM calculation the recently recalcu- 
lated cross sections [13] ■ 

small Doppler shifts. This reveals that cross sections re- 
sponsible for the development of such low energies coming 
from an interplay of Coulomb de-excitation and elastic 
collisions could still be deficient. In addition, the cascade 
calculation starting at n — 8 possibly neglects stronger 
effects from Coulomb de-excitation in the outer shells 
(n > 11), where the ESCM uses a classical trajectory 
Monte-Carlo approach. 

Consequently, as a further approach a "model free" 
method was applied, which has already been used in the 
case of the neutron TOF spectra [8j . The Doppler contri- 
butions from Coulomb de-excitation to the /xH(3p — Is) 
line shape were determined by modelling the kinetic 
energy distribution with several rectangular boxes, the 
number of which is assessed by the transitions in the 
preceding cascade. 

Several sets of kinetic energy boxes were investigated 
including besides An = 1 also An = 2 Coulomb tran- 
sitions. Hyperfine splitting, i.e., both line positions, 
relative population, (flat) background, and the relative 
weight of the boxes were free parameters in the fit. The 
sum of the relative weights of the boxes was always nor- 
malized to one. 

It is a major result of this study, that three nar- 
row boxes are essential but also sufficient to model a 
kinetic energy distribution yielding a good fit to the 



/iH(3p — Is) line shape: (i) one low-energy component be- 
low Tkin = 2 eV collecting fiH atoms which gained their 
energy from n > 10 Coulomb transitions and/or high 
velocity systems degraded by collisions, and (ii) two 
at higher energies corresponding to the de-excitation 
steps n = 5 -> 4 (T kin = 26.9 eV) and n = 4 -> 3 
(T fem = 58.2eV). 

A x 2 analysis shows that a weight of « 60% is manda- 
tory for the low-energy component. In case of the high- 
energy components the fit is only sensitive to the rela- 
tive weight and the central value. Extending the bound- 
aries up to ± 30% of the central values, affects the re- 
sult by less than 1.4 standard deviations. Therefore, the 
kinetic energy distribution could be condensed to three 
narrow intervals. 

The best reduced \ 2 is found to be \r = 0.947 for 
the kinetic energy intervals set to [0-1.8], [26.4-27.4], and 
[57.7-58.7] eV resulting in relative weights of (61±2)%, 
(25±3)%, and (14±4)% (FigM "best fit" in the "model 
free" approach). Uncertainties represent statistical errors 
only. 

A correlation study of hyperfine splitting and relative 
population was performed by using the three kinetic en- 
ergy intervals found in the above mentioned analysis, but 
with their weights, total intensity, and background kept 
as free parameters. The best Xr — 0.941 is obtained 
for a hyperfine splitting of 211 ± 19 meV, a triplet to sin- 
glet population of (3.59 ± 0.51) : 1 (Fig.[3]-A), and relative 
weight of 61 ± 1% for the low-energy component, where 
errors correspond to la. The 1, 1.5, and 2 a contours are 
also shown. When fixing the hyperfine splitting to the 
theoretical value, the best fit is obtained for a triplet 
to singlet population of (2.90 ±0.21) : 1 (Fig.[3j-B), very 
close to the statistical value. The x 2 differs only by 1.5 a 
from the best value. 
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FIG. 3: x 2 contour for the correlation of hyperfine splitting 
and relative population in the "model free" approach. The 
dashed-dotted line displays the location of the minimum x 2 
for the corresponding hyperfine splitting. (A - D: see text). 
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The results from the present experiment led to a recon- 
sideration of cross sections involved in exotic-hydrogen 
de-excitation. In a fully quantum-mechanical close- 
coupling approach, elastic scattering, Stark transitions 
and Coulomb de-excitation now have been calculated in 
a unified manner 0, @, [26| . 

Including these cross-section results in the ESCM 



model [13] , the relative weight of the low energies (T> 



kin. 



< 



2eV) increases to 55%, which is closer to the constraint 
found in the "model free" approach. A fit to the mea- 
sured /iH(3p — Is) line shape by using directly the new 
kinetic energy distribution (Fig. [J) yields a much bet- 
ter description of the data (Fig. [2]). Leaving hyperfine 
splitting and triplet to singlet intensity ratio as free pa- 
rameters, values were obtained of 194 ± 12meV for the 
hyperfine splitting and (2.94 ±0.24) : 1 for the triplet to 
singlet intensity ratio (Fig.[3]-C). The striking agreement 
is indicated by \r — 0.997. When fixing the splitting to 
the theoretical value the relative triplet to singlet popu- 
lation becomes (3.17± 0.27) : 1 (Fig.dfD). 

Another process, which may affect the line shape, is 
due to molecule formation with subsequent radiative de- 
oay (np) n i + H 2 -> [{^pp)*pee\* -> [(fipp)p ee]* + 7 
27J,|28|. If radiative decay contributes significantly, a line 
broadening or even satellites at the low-energy side are 
expected. In the case of [ipp, a branching ratio for radia- 
tive decay was predicted at the few per cent level 29|, 3(| • 
However, in this experiment, no evidence was found at 
the 1% level for any broadening except from Coulomb de- 
excitation. The result is corroborated from the absense 
of any density dependence of the 7rH(3p — Is) transition 
energy 0]. 

To summarize, the line shape of the /TI(3p — Is) tran- 
sition was measured with a high resolution Bragg spec- 
trometer. The influence of Coulomb de-excitation was di- 
rectly seen from a line broadening compared to the spec- 
trometer resolution. By using a "model free" approach 
various Doppler contributions are identified, which are 
attributed to preceding Coulomb transitions. A large 
fraction of the /iH systems are found to have kinetic ener- 
gies below 2 eV. The measurement yields the /iH ground 
state hyperfine splitting as calculated from QED and 
confirms experimentally the statistical population of the 
triplet and singlet Is states. The measurement triggered 
a new calculation of cross sections resulting in a signifi- 
cantly improved description of the /iH(3p— Is) line shape 
and serves as a basis for a further evaluation of the de- 
termination of the isovector scattering amplitude from 
pionic hydrogen data with increased precision. 
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